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ABSTRACT 

We perform time-resolved spectroscopy of all the type I bursts from the Rapid Burster 
(MXB 1730—335) detected with the Rossi X-ray Timing Explorer. Type I bursts are 
detected at high accretion rates, up to ~ 45% of the Eddington luminosity. We find 
evidence that bursts lacking the canonical cooling in their time-resolved spectra are, 
none the less, thermonuclear in nature. The type I bursting rate keeps increasing with 
the persistent luminosity, well above the threshold at which it is known to abruptly 
drop in other bursting low-mass X-ray binaries. The only other known source in which 
the bursting rate keeps increasing over such a large range of mass accretion rates is 
the 11 Hz pulsar IGR J17480— 2446. This may indicate a similarly slow spin for the 
neutron star in the Rapid Burster. 

Key words: stars: neutron - X-rays: binaries - X-rays: bursts - X-rays: individual: 
MXB 1730-335 



1 INTRODUCTION 



currence time, energy and spectral properties (for a review. 



First discovered by Lewin et al. (19761, the Rapid Burster 



(MXB 1730-335, hereafter RB) is a neutron star (NS) low- 
mass X-ray binary (LMXB) located in the globular cluster 
Liller 1, with a recent distance measurement of 7.9 ±0.9 kpc 
( Valenti, Ferraro fc Origlia|2010[ ). It is a recurrent transient, 
which had very regular outbursts every ~ 210 d for most 
of the time up to MJD ~ 51500, after which the outburst 
recurrence time abruptly dropped to ~ 100 d ( ||Masetti^2002'| ). 

The RB is unique, in that it features both type I and 
type II bursts. While the former are due to the heating and 
cooling of the surface after a thermonuclear flash in the un- 
derlying layers (for a recent review see [Galloway et al.|2008| 



henceforth G08), type II bursts are most likely due to the 
release of gravitational energy from the inner accretion disc 



during sudden accretion events (as first pointed out by Hoff- 
man, Marshall fc Lewin|1978 l, probably linked to magnetic 
gating of the accretion flow (Lamb et al. 1977 Spruit & 



Taam||1993| [D'Angelo fc Spruit||2010| ). Consequently, these 
two categories of bursts can differ greatly in duration, re- 
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Lewin, van Paradijs fc Taam|1993 hereafter LPT93). 
Typically, type I burst recurrence times tree are of the 
order of hours. Although some type I bursts have been re- 
ported to recur within ~ 200 s (e.g. [Keek et al.||2010[ l only 
one source - IGR J17480-2446, hereafter T5X2 - exhibited 
that behavior for prolonged series of bursts lasting for hours, 
instead of multiple events consisting just of two, three or 
four bursts. [Linares et al.| ( |2012[ ) (henceforth L12) reported 
T5X2 to constantly increase in burst rate as its persistent 
luminosity rose from 0.1 to 0.5 of the Eddington luminos- 
ity I/Edd, with bursts transitioning to a mHz quasi-periodic 
oscillation (QPO), with a period as short as 240 s. 

This behaviour is well in agreement with models for 
thermonuclear bursts. Theory predicts the bursting rate to 
keep increasing with the persistent luminosity almost up to 
the Eddington limit, w hen stable burning sets in ([Fujimoto,[ 



Hanawa fc Miyaji|19"8T||Heger, Gumming fc W oosley |2007 1 
still, the burst rate for the majority of burst sources is ob- 
served to decrease dramatically above a few per cent of the 
Eddington luminosity (e.g. Gornelisse et al. 2003 Fig. 16 
in G08). Although the basic physics of type I bursts is well 
understood, this discrepancy between the expected and ob- 
served troc remains an unsolved issue in our understanding 
of the phenomenon. 

The unusual properties of T5X2 probably set it apart 
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from the other bursting sources, and could possibly explain 
its extremely regular bursting behavior (L12). T5X2 is an 
11 Hz pulsa r with a reported magnetic field betwe en 10* 



and 10^° G ( jPapitto et al.||201l[ [Miller et al.||201l| ), mak- 



ing it the slowest rotating burster known, and possibly one 
of the high magnetic field bursters as well. Slow rotation 
could prevent turbulent mixing from stabilizing burning at 
low mass-accretion rates ( Keek et al.[ 2009') and reduces the 
influence of the Coriolis force on flame spread and confine- 
ment ( Spitkovsky et al.|2002 I. Magnetic confinement, which 
probably accounts for the presence of burst oscillations in 
T5X2 ( Cavecchi et al.||2011 l, could limit burning to a con- 
stant portion of the stellar surface, while the burning area 



may vary with mass-accretion rate in other sources ( Bildsten 



2000) 



As already mentioned, the presence of type II bursts 
is probably indicative of a prominent dynamic role of the 
magnetic field B in the RB. We therefore seek to investi- 
gate whether its type I burst behaviour shows analogies to 
T5X2 that could be explained in terms of a slow spin or of 
surface magnetic field effects, related to a magnetic gating 
mechanism in the accretion flow. 

We perform a systematic study of all the type I bursts 
from the Rapid Burster detected with the Proportional 
counter array (PCA) onboard the Rosst X-ray Timing Ex- 
plorer (RXTE), increasing the number of bursts over that 
already known by a factor of two. We observe bursts over one 
of the largest ranges in persistent flux, and interestingly find 
that the burst rate increases to high values with the mass 
accretion rate. We discuss the implications of this behaviour 
in our understanding of the peculiar RB. 

Section [2] describes the data selection, with particular 
care devoted to how we excluded spurious bursts (type II, or 
type I from a different source). Section [s] explains the anal- 
yses, spectral and morphological, performed on the bursts 
and the persistent emission. In Section |4] we report our re- 
sults on the burst behaviour as a function of the persistent 
luminosity of the RB. Finally, we discuss the implications 
of our findings for the burst regimes and the source spin in 
Section (5] 



2 OBSERVATIONS 



2.1 Dataset 



The PCA ( [Jahoda et al.|[2006| ) consists of five co-aligned 
proportional counter units (PCUs) that combine to a total 
effective area of 6000 cm^ at 6 keV, in a 2 to 60 keV band- 
pass. The photon energy resolution is 18% full width at half 
maximum at 6 keV and the time resolution is programmable 
down to 1 fis. The field of view (FOV) is circular, with radius 
~ 1° (full width to zero response). 

The RB has been extensively monitored with the PCA, 
from the beginning of the mission in early 1996 until its 
shutdown in January 2012. We collected all PCA data that 
we could find when the RB was in the FOV of an active 
PCU, for a total exposure time of 2.83 Msec (this is the sum 
of all so-called good-time intervals, excluding times when 
RXTE is slewing, near to an SAA passage, or experiences a 
large particle rate) . 

A major issue of the PCA data is that the RB has an 



angular separation of only 0.56° to the persistently active 
burster 4U 1728-34. In fact, most observations (67%) were 
aimed at 4U 1728-34. The RB was the aimpoint for 14% 
of the time, so that both sources are in the FOV for 82% 
of the time (including as well 1% when neither source was 
at the centre of the FOV). During the remaining 18%, the 
PCA aimpoint was offset so that only the RB is inside the 
FOV, decreasing the collimator response by a factor ~ 0.4, 
but allowing for spectroscopic measurements unaffected by 
contributions from other sources. We note that, due to the 
transient nature of the source, only 32% of the data (or 
0.91 Msec) was taken whilst the source was actually active. 

2.2 Identification of RB type I bursts 

We employed the Standard 1 data product, consisting of 
0.125-s resolved photon count rates per PCU without any 
photon energy resolution, to search for bursts. Our method 
consists primarily of a computer algorithm to search for 
statistically significant upward fluctuations above a steady 
background on time-scales of 1 to 300 s. We verified the re- 
sult of this search by a careful visual inspection of the light 
curves and the identified bursts. We additionally applied the 
criterion that the Rapid Burster should be active, as may be 
verified with independent RXTE-ASM or ISS-MAXI mea- 
surements, to exhibit bursts. The computer algorithm was 
set to accept fluctuations above the background if they have 
a significance accumulated over the burst that is at least 10 a 
and if they are accompanied by a stable pointing of the PCA. 
The algorithm is efficient as long as there are time intervals 
before or after a burst within 300 s where the background 
can be measured confidently. As a result some long type II 
bursts may be missed that have no flat background between 
them and bursts that occur during slews with the RB in the 
field of view. Through a visual inspection, however, we made 
sure that type I bursts with the typical fast-rise exponential- 
decay profile were included even if they occurred during a 
slew. We find a total of 7261 bursts, including 121 type I 
bursts and 7140 type II bursts. The latter number is a lower 
limit. 

Type I bursts from 4U 1728—34 are easily recogniz- 
able by their short durations (< 10 s, typical of H-poor 
nuclear fuel) and by their brighter (often Eddington lim- 
ited) peak fiuxes, bimodally distributed about a mean of 
9.2xl0~* erg cm~^ s~^ for the bursts showing photospheric 
radius expansion (PRE) and 4.5x 10~® erg cm~'^ s~^ for nor- 
mal bursts ( |Fox et al.|[2001[ [Galloway et al.||2003| G08), a 
factor of 6 brighter than most bursts from the Rapid Burster. 
For an angular separation of 0.56°, this ratio decreases by 
40% when the Rapid Burster is at the centre of the FOV and 
there are some type II bursts from the Rapid Burster which 
come close to half the peak flux of bursts from 4U 1728-34. 
Still, these are rare and the bursts can be distinguished by 
their time proflle. 

It is also straightforward to exclude most type II bursts. 
Long-lasting, fiat-topped bursts, preceded and followed by 
dips in the persistent emission are clearly recognizable as 
type II, as are the intermediate-duration bursts showing 
multi-peaked decays (Tan et al. 19911. The shortest type 



II bursts, the only ones featuring single, short peaks and a 
fast rise, are potentially ambiguous. They however have very 
short (< 10 s) decays and recur very fast in series of 8-40 
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bursts, making their cumulative energy release incompatible 
with a thermonuclear origin (LPT93, see Secjs] for further 
discussion) . When a type I burst hides in a sequence of such 
bursts it is easily recognized by its much slower decay. 

We nevertheless acknowledge that we could have erro- 
neously misidentified a few type I bursts as type II. We in- 
tend to publish a comprehensive overview of all type I and 
II bursting behavior of the Rapid Burster elsewhere, with 
a detailed description of the burst identifications. Here we 
concentrate on the 121 type I X-ray bursts that we identi- 
fied from the Rapid Burster in the PCA data. 65 have al- 
ready been reported ( Fox et al.|2001 G08). Of these bursts, 
some have only been partially observed, due to their prox- 
imity to the beginning or end of a pointing, and some are 
surrounded or even covered by trains of fast-recurring type 
II bursts. These bursts are not suitable for analysis, espe- 
cially when the persistent emission is very low, because part 
of the accretion is then taking place through the type II 
bursts, implying that in this case the persistent flux does 
not even qualitatively reflect the actual mass-accretion rate 
on the surface, and a relation between burst properties and 
the persistent emission should therefore no longer be ex- 
pected. Consequently, we only focus on type I bursts from 
observations in which no type II bursts were observed. 



3 ANALYSIS 

We set out to determine a number of key parameters for all 
bursts. 



3.1 Burst durations and recurrence times 

Using the 1-s resolution light curve, we modelled the tail of 
the burst light curve, when the intensity is less than 90% of 
the peak intensity minus the pre-burst level, with one or - 
when necessary - two exponential decay functions. We used 
these fits to determine the burst duration idur, defined as the 
time difference between the burst start time (see G08) and 
the time when the intensity drops below 10% of the peak 
intensity minus the pre-burst level, when the burst count 
rates are still significantly above the noise, so that variations 
in the persistent count rate do not affect our estimate of 
the burst duration. Generally speaking, single exponentials 
provide sufficiently good fits to shorter bursts (tdur ~ 50 
s), while two exponentials are necessary for longer bursts 
(tdur ~ 100 s). Although such a definition of the burst du- 
ration is not affected by the varying noise level at different 
persistent count rates, the burst emission is still not over at 
tdur. Spectra were therefore extracted up to the time tend 
when the intensity decreases to within 2o of the pre-burst 
level. We then used the exponential fits to estimate how 
much fiuence of the burst may have been contained in the 
far tail of the light curve (i.e., after tend). This factor, be- 
low 5% for most bursts, was added to the burst fiuence and 



to its uncertainty (see Sec 3.2.21, to avoid a systematic ef- 



fect towards less energetic bursts at higher persistent count 
rates. 

The burst recurrence time tree is defined as the time 
elapsed between the start time of the previous burst and that 
of a given burst. When data gaps occurred, the recurrence 
time was labelled an upper limit. In five cases, however. 
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Figure 1. Cooling type I bursts. Notice the slightly different x 
and j/-axis ranges. The bursts are very regular, with separations 
varying by just 10% in the upper sequence (as the average -Fpers 
decreases by roughly the same amount) and 5% for the lower se- 
quence. The burst separations covering the data gaps are close to 
twice the mean separation of the other bursts, suggesting another 
burst occurred in the gap (see Sec. |3.1| . In the lower panel, the 
last burst took place in a configuration that excluded 4U 1728-34 
from the FOV and reduced the RB intensity by a factor ~ 0.4. 
We rescaled the count rate accordingly. The remaining difference 
in the persistent emission level is therefore due to the lack of 
contribution of 4U 1728-34. 



the observed tree was very close to being an integer multiple 
(~2 or 3 times longer) of the recurrence time observed in the 
adjacent bursts, which was seen to be stable over series of up 
to five bursts, and a burst occurring with the average value 
of tree would have fallen in a data gap (see Fig.[T]). Therefore, 
we assumed these bursts to be recurring on a regular time- 
scale and divided their observed tree by the integer amount 
of times that gave a nearly constant burst rate. 



3.2 Burst energetics 

Throughout our analysis, the particle and cosmic back- 
ground, as determined with ftool pcabackest, were sub- 
tracted from all spectra. Response matrices were generated 
with pcarsp (v 11.7.1). All active PCUs were employed, and 
a correction was applied for deadtime (although it is always 
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small for the RB). Low-energy absorption by the interstellar 
medium was taken into account using the model of |Morrison| 
& McCammon (19831, with an e quivalent hydrogen column 
density of iVe = 1.6 x 10^^ cm"^ (iFrogel, Kuchinski & Tiede 



1995[|Masetti et al.|2000l|Falariga et al.|2004[ ). Following the 
RXTE Cookbook prescriptiorlM the analysis was limited to 
the 3 — 30 keV range, and a systematic error of 0.5% was 
added to the results. 

To determine the energetics of the bursts and of the per- 
sistent emission, we performed a time-resolved spectroscopic 
analysis according to the following procedure. 



3.2.1 The persistent emission 

First, the spectrum of the non-burst X-ray emission was ex- 
tracted from 496 s long Standard 2 data stretches preced- 
ing the burst, to constrain the background signal underlying 
the burst. We modelled this spectrum with a generic LMXB 
model consisting of a disc black body ( Mitsuda et al.|1984 |, 
a power law and a Gaussian centred at 6.4 keV (G08), that 
yielded a Xred very close to one for all the pre-burst spectra. 
In the many cases that the FOV also included 4U 1728—34, 
the measured flux also included a contribution from that 
source and the resulting fit parameters are not a truthful 
representation of the persistent emission of the RB, but for 
assessing the burst energetics this is not crucial as long as 
the model is a good description of the non-burst emission. 

For those cases where 4U 1728-34 was outside the FOV, 
this model also provides a good measurement of the per- 



tamination was an issue, a dedicated offset observation (see 
Sec. 2.1 1 taken no longer than one hour before or after the 



burst was often available. We used these uncontaminated 
spectra to estimate the persistent flux Fpers in the range 3- 
25 keV, resorting to the same three-component model men- 
tioned earlier. All but six such measurements were taken 
within 1500 s from the burst onset. In all cases, we took 
care to verify the stability of the average count rate over the 
interval between the measurement and the burst. The count 
rate variation in the RB over the course of an hour was al- 
ways below 5%. All fits were satisfactory, giving Xred ^ 1-5 
for 40-50 degrees of freedom. Note that in the case of short 
recurrence time bursts, a single -Fpors measurement may ap- 
ply to multiple bursts. 

For each of these uncontaminated observations, a bolo- 
metric correction Cboi was calculated from the portion of 
the observation excluding 4U 1728—34. We fitted an ab- 
sorbed Comptonisation model (compTT in XSpec) to the 
spectrum extracted for each PCU in the energy range 3- 
25 keV, and integrated the best fit model to give the flux 
over this range. We then created (in XSpec) a "dummy" 
(ideal) response covering the broader energy range of 0.1- 
200 keV, and calculated the integrated, unabsorbed flux over 
that range. The bolometric correction for a given PCU was 
estimated as the ratio of these two fiuxes, and the value 
adopted for each observation was the mean of the values for 
the active PCUs, with the standard deviation of the values 



^ See |http : //heasarc ■ nasa ■ gov/docs/xte/reclpes/cook_book ■ ] 
[EEiI| 



over all PCUs adopted as the uncertainty. The Comptonisa- 
tion model was chosen over the phenomenological (disk black 
body plus power law) because the Comptonisation model 
has a finite integral for any choice of the input parameters. 
In fact, because the persistent spectrum was generally soft, 
the electron scattering temperature kT^ was typically a few 
keV, and could be well constrained by the spectrum up to 
25 keV. For a few of the observations, the spectrum was 
harder, and in some cases the kTe could not be constrained 
by the fit. For most of those observations, the addition of 
a Gaussian component improved the fit sufficiently that the 
kTe could be constrained; for a few examples it was possible 
to constrain the plasma temperature for all but PCU #3, 
and for those observations that spectrum was excluded. 

We adopt FpcTs = CboiJ3-25 kov as the bolometric inten- 
sity for the source in each observation. We also often refer 
to its ratio to the Eddington flux FEdd. To derive the latter, 
we apply the distance estimate to the Eddington luminosity 



LEdd = 3.5 X 10^ 



M 



1.4M, 



1 + X 



ergs s 



(1) 



with M the stellar mass and X the hydrogen fraction (e.g. 
G08). This corresponds to an Eddington flux from the disc 
i^Edd = 2.8xl0~® erg cm~^ s~^ for a "canonical" NS with 
mass M = 1.4Mq and radius _R = 10 km, accreting solar- 
composition material {X — 0.7). 

3.2.2 The bursts 

The burst data were divided into a number of time bins, 
varying in duration to keep the photon count, and hence 
the relative error on the derived quantities, approximately 
constant. The bin duration was constrained to be at most 1 s 
during the burst rise, to be able to track the initially rapidly- 
varying temperature. The minimum bin duration possible 
was constrained by the time resolution of Standard 1 data, 
at 0.125 s. 

Spectra were extracted for each time bin from event 
files, the time resolution of which varied between 1 and 
125 /is. The non-burst emission was not subtracted. Spec- 
tra during the bursts were modelled by the combination of 
black body radiation (leaving free temperature and emis- 
sion area, to account for the varying burst emission, and 
multiplied by a model for interstellar absorption with TVh ~ 
1.6 X 10^^ cm~'^) and the model as found for the non-burst 
spectrum discussed above, keeping the parameters of the 
latter fixed to the pre-burst values. 

The assumption underlying this "standard" procedure 
is that Fpeis is constant during the burst. The well-known 



caveat (van Paradijs & Lewin 19861 is the possibility that 



the persistent emission actually varies during the burst, giv- 
ing rise, for instance, to systematic errors for the radii that 
can get very large in the burst tails. Variations could arise 
because of radiation effects on the inner accretion disc, or 
because of the NS becoming sufficiently hot to contribute 
to the persistent spectrum. Worpel, Galloway & Price (sub- 
mitted) noted a rise in persistent emission during radius- 
expansion bursts which they attribute to radiation drag ef- 
fects. We performed a similar analysis, adding a multiplica- 
tive factor /a by which the persistent emission is multiplied, 
as a variable parameter in our spectral fits. However, even 
for the most luminous burst, the effect appeared negligible. 
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with /a — 1.06 ± 0.29 at the burst peak. It therefore seems 
unhkely that this is a major factor for RB bursts, which are 
all significantly sub-Eddington. 

The black body description of our burst spectra is al- 
ways satisfactory. This simple, physically founded model 
provides a fairly accurate means to obtain bolometric fluxes 
and fluences as follows. The black body fit to each mea- 
sured spectrum results in the measurement of two parame- 
ters: the temperature, expressed in kT, and normalization 
N, expressed in terms of the square of the radius in km of 
a sphere assumed to emit the black body radiation at a dis- 
tance of 10 kpc. The values of both these parameters come 
with a 1-sigma uncertainty which defines a 68% confidence 
margin for a one parameter fit. The burst bolometric fiux is 
calculated from Stefan-Boltzmann law: 



Fbb = 1.07 X 10" 



TV 



kT Y 



(2) 



To determine the uncertainty in Fhh, we sample 10,000 
values of kT and N from Gaussian distributions centered at 
the fitted values and with standard deviations equal to the 
respective 1-sigma uncertainties, calculate for each sample 
-Fbb and from the 10,000 samples calculate the mean and 
standard deviation. The fiuence is calculated from the sum 
of all bolometric fluxes times the integration times of the 
related spectra and the uncertainty in the fiuence from the 
root of the quadratic sum of the uncertainty in all terms. 
Since the errors in kT and A'^ were treated as if these pa- 
rameters are independent, the error in the burst bolometric 
flux can be considered as conservative. 

Some uncertainty remains for the fiuence, for the burst 
radiation in the possibly long tail where it drowns in the 
noise. We estimated this tail contribution from the expo- 
nential fits discussed in Sec. 13.11 



3.2.3 Determination of a and (3 



The ratio of the integrated bolometric persistent fiux to the 
burst fiuence a = Fpcrsircc/-Bb could be precisely determined 
only for sixteen bursts. The other bursts suffered from source 
confusion (see Sec. |2.1[ | or data gaps making tree only an 
upper limit. For five additional bursts, however, we could 
establish the recurrence time fairly confidently even in the 
presence of data gaps (see Sec. 3.1 1, bringing the total num- 



ber of bursts with a known value for a to 21, including four 
bursts lacking spectral evidence for cooling in their decays 
(see Sec. 4.1 1. Errors on Fpcrs and Sb are nearly always be- 



low 10%, and the resulting errors on a are generally below 
20%. 

We also looked at the ratio of the burst peak to per- 
sistent fiux /3 = -Fpcak/Jpors, which we could derive for 57 
bursts. The errors on the peak flux were somewhat large, 
which propagated to errors on /3 getting large especially at 
the lower end of the fpcrs range. 



4 RESULTS 

4.1 Lack of cooling 

We fail to find evidence for cooling in six bursts, meaning 
the temperature difference AT = Tpcak — Tkur is consistent 



with zero within error bars, with the black body tempera- 
ture Tpoak measured at the ffux peak and Tdur at the time 
tdur. They occurred at the beginning of the June- July 1997 
outburst, at the highest level of -Fpcrs we observed in all 
data (see Fig. [2| . The highest value for which cooling is not 
observed is /3 = 0.37 ± 0.09, while the lowest value among 
bursts with a cooling tail is /3 = 0.50 ± 0.10. On average 
P = 0.30 ± 0.13 for the non-cooling bursts. Their a values 
range from 87 ± 13 to 153 ±26, consistent with a thermonu- 
clear origin, and they have time profiles typical for type I 



bursts (see Sec 2.2 1. This convinces us that they are not type 
II bursts, for which reported values of a range from 0.15 to 
4.29 ( |Tan et al.|199H|Masetti et al.|200"o| . 



4.2 Burst behaviour as a function of accretion 
rate 

In Fig. |3] all burst parameters are plotted against the bolo- 
metric accretion fiux Fpava and its ratio to the Edding- 
ton fiux ^Edd- The RB produces regular type-I bursts over 
a much wider range of persistent flux than seen in other 
systems, from (2.7 ± 0.2) x 10"^ erg cm'^ s'^ to (12.5 ± 
0.3) X 10"^ erg cm"^ s"\ or between 10% and 45% of F-em- 
Some trends in burst parameters versus persistent flux 
are clearly visible in Fig. [3] The ratio (3 clearly decreases with 
increasing Fpers, until it falls below the threshold /3 = 0.7 



which Linares, Chakrabarty & van der Klis (20111 suggest 



determines whether cooling can be detected in PCA data 



(Linares, Chakrabarty & van der Klis 20111. This is not 



only trivially due to the increeise in Fpers- The peak fluxes 
show a decreasing trend with increasing Fp^is, down from 
a maximum of fpeak = (13.6 ± 2.6)xl0~^ erg cm~^ . 
On the other hand, the non-cooling bursts at the high- 
est -Fpers show much lower peak fluxes, averaging Fpeak — 
3.7x10"^ erg cm"^ s"\ 

Most importantly, the recurrence time tree falls steadily 
with increasing Fpcrs, down to a minimum of 467 s. Re- 
gardless of -Fpors, the bursts appear to be very regularly 
clocked (see Figs, [l] and |2|. A fit to the recurrence times 
yields a power law with index —1.81 ± 0.04. The fit how- 
ever is unacceptably poor, with x?od — 22.1 for 22 degrees 
of freedom. Excluding the non-cooling bursts, this becomes 
a fiatter —0.95 ± 0.03, with x?ed = 3.46 for 17 degrees of 
freedom. The impossibility of fitting a single relation is due 
to the steepening of the relation in the non-cooling regime, 
as is clear from Fig. [3] However, too few non-cooling bursts 
are available and they span too narrow an Fpcrs range for a 
fit to this latter regime to be meaningful. 

Despite the large spread in values, both iJb and tdm 
seem to stay relatively constant until the highest persistent 
fluxes. A drop then occurs, coincidently with the disappear- 
ance of cooling from the burst spectra. One might suspect 
this to be due to the poorer statistics in the burst tail at 
large -Fpcrs - a sizeable fraction of the burst fiuence might be 
lost in the noise. However we defined tdur and E^, exactly to 



avoid such an effect (see Sec. 3.1 1. We therefore believe this 



drop to be genuine. The average rise time shows no corre- 
lation with the persistent fiux (triso — 7.1 s, with standard 
deviation a ~ 2.1 s), although its average value is different 
for the cooling and non-cooling bursts (6.8 and 9.6 s respec- 
tively), perhaps hinting at a late change in bursting regime. 
Together with the change in burst duration tdur, which on 
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Figure 2. Fast recurring, non-cooling type I bursts. Concerning 
the oflFset in count rate around the last burst, see Fig. [T] 



the other hand seems to get shorter at the largest accre- 
tion rates, this makes for bursts becoming somewhat more 
symmetric in shape. 

Excluding the bursts at the highest accretion rates, 
nearly constant a values are observed for the cooling bursts, 
although at lower Fpcrs the longer tree means that often only 
upper limits on a can be measured. For the cooling bursts 
for which a is well constrained, on average a ~ 46 with a 
standard deviation of ~ 10. Including the bursts for which 
we divided tree by an integer value estimated by comparison 
with the surrounding bursts (see Sec. |3.1| and Fig. [T| yields 
an average a ~ 44 with a standard deviation of ~ 9, reinforc- 
ing our inference that our assumption of stable recurrence 
time is correct. Only at the high-end of the distribution is 
an upturn in ct visible, peaking to a = 153 ± 26. On av- 
erage, the non-cooling bursts have a ~ 98 with a standard 
deviation of ~ 26. It is important to stress that for all the 
bursts for which a is well constrained, the value points to a 
thermonuclear origin. 

No obvious explanation is available for the outlier, the 
very low Fpeak burst at Fpers = 3.6±0.3x 10~^ erg cm~^ s~^. 



5 DISCUSSION AND CONCLUSIONS 

5.1 Cooling and burst type identification 

The identification of a burst from the RB as type I or II 
can be difficult. Several attempts have been made using 
broad-band spectroscopy, but the results are inconclusive 



at best (Kunieda et al. 


1984 Stella et al. 


1988 


|1991| |Lubin et al.||1992 


|Guerriero et al.| 


19991) 



Tan et al. 



less, the presence of cooling in time-resolved spectroscopy 
of type I bursts has long been accepted as a hallmark of 
thermonuclear burning (LPT93), for type II bursts generally 
exhibit very little spectral evolution. This has been used to 
argue for the occurrence of type II bursts in other sources 
(e.g. T5 X2; IGalloway fc in 't Zand||2010[). However several 
authors ( Motta et al.||20lT Chakraborty & Bhattacharyya 



2011) soon pointed out that the a values for T5X2 bursts 



that did not show a cooling trend were consistent with a 
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thermonuclear origin. Linares et al. (2011) found the occa- 



Figure 3. Peak flux -Fpcaki ratio /3, fluence E\^, burst duration 
tdur, burst recurrence time tree and ratio a plotted against the 
bolometric unabsorbed -Fpors (lower x-axis) and its ratio to the 
Eddington flux (upper x-axis). Horizontal error bars refer to the 
lower x-axis. Empty and filled symbols correspond to cooling and 
non-cooling bursts respectively, triangles indicate upper or lower 
limits, squares refer to tree values (and the corresponding a val- 
ues) that were inferred dividing the observed value by an integer 
to match the ones in the nearby bursts (see Sec. |3.1| . For the 
"cooling threshold" at /3 = 0.7, see Sec. |4.2| The fit to tree is for 
cooling bursts only. 



Bursts from the RB: indications for a slow rotator 7 



sional lack of cooling in some bursts from T5X2, Cir X-1, GX 
17+2 and Cyg X-2 to coincide with a low value for /3 near the 
outburst peak. They proved the a values of these bursts to 
be consistent with a thermonuclear origin, and found cool- 
ing to disappear from bursts below a ratio /3 of 0.7. They 
attributed the lack of measurable changes in Thh to two fac- 
tors: a loss of sensitivity due to the increased luminosity of 
the persistent emission (resulting in lower signal to noise); 
and a reduction in the column depth required for ignition 
as the layer temperature rises (resulting in less mass being 
burned and the bursts reaching lower peak temperatures). 
Values of f5 for non-cooling PC A bursts from Cir X-1, GX 
17-1-2 and Cyg X-2 are also below /3 ~ 0.7, ruling out these 
sources as type II burst emitters. 

Ultimately, a remains the most reliable observable to 
distinguish type I from type II bursts. Depending on the fuel 
composition, thermonuclear burning can release between 1.6 



and 4.4 MeV/nucleon (Wallace & Woosley 1981 Fujimoto 
et al.|1987 l, while the energy release from accretion on to a 
"canonical" neutron star is about 180 MeV/nucleon, yield- 
ing expected a values between ~ 40 and ~ 110 for thermonu- 
clear bursts, about one to two orders of magnitude above 
what is observed for type II bursts ( Tan et al.|1991 Masetti 
et al.|[2000 l. We found compelling evidence that RB bursts 
at the highest persistent fluxes have a thermonuclear origin 
with precise measurements of their a values, and showed 
that cooling disappears below a /3 threshold that is roughly 
consistent with that inferred by Linares et al. (20111. Al- 



though the transition seems to lie at a somewhat lower value 
of /3 in the RB, it should be noted that the energy band they 
used to calculate Fpers is narrower than ours (2 — 50 keV ver- 



sus 1 — 200 keV). Also, Guerriero et al. ( 1999 1 already noted 



that in the RB type II bursts only appear late in the outburst 
decay, usually after the persistent luminosity has decreased 



to Fp, 



3x10 erg cm s . For all but one burst we 



analysed, Fpors lies above this value and we confirm that no 
type II bursts are observable in the lightcurves. It is there- 
fore unlikely that type II bursting activity resumes at the 
highest -Fpcrs levels, with a values two orders of magnitude 
larger than usual and smooth exponential decays rather than 
the multi-peaked decays type II bursts of similar duration 
are reported to show ( |Tan et al.|1991[ ). 

The paucity of identified non-cooling bursts is due not 
only to the relatively short time the source spends at the 
highest accretion rates during an outburst, but also to the 



drop in outburst peak flux after MJD 51500. Masetti (20021 



reported a drop in both the outburst peak flux and recur- 
rence time by roughly by a factor of two after that date. 

5.2 Bursting regimes 



We showed that the RB keeps emitting type I bursts at 
an ever increasing rate, over a range of -Fpors from (2.7 ± 
0.2)xl0"^ erg cm"^ s"^ to (12.5±0.3) x 10"^ erg cm~^ s"^ 
Assuming a distance of 7.9 kpc, this corresponds to (2.6 — 
9.3)xl0^^ erg s-\ or (10 - 45)% LEdd- The X-ray flux is 
howe ver known to b e a poor proxy of the real mass accretion 
rate d van der Klis||2001) , so that the real range in M that 



the source is spanning remains uncertain. Accounting for the 
uncertainty in the distance, the maximum luminosity at the 
time of a burst ranges between 35% and 53% of Lfidd- 

That the bursting rate keeps increasing over such a large 



range in Fpcrs is a striking similarity between the RB and 
T5X2, and it sets these two sources apart from all other 
known bursters (e.g. Cornelisse et al.|2003 Fig. 16 in G08). 
The trend in Fpeak is also similar. As mentioned, we do not 
observe a continuous decrease in Eh and tdur- These trends 
are, however, rather dubious in T5X2 when looking at the 
individual bursts, and only become clear switching to the 
daily averages of these burst properties (see Fig. 7 and 9 in 
L12). A larger number of available bursts in T5X2 is due 
to a bursting rate that initially grows faster: Ucc oc -Fp^^s 
0.3, after which it settles on a flatter 



up to Fpers /-Ffidd 

tree OC -Fj^rs. The latter relation seems closer to the one fol- 
lowed by the RB for all but the highest Fpers levels: exclud- 



ing non-cooling bursts, we fit tree oc Fj 



(-0.95±0.03) 



pcrs 



This is 



consistent with the empirical relation s tree oc Fp^rs found 
in GS 1826-24 (ICalloway et al.||2004|| and troc oc Fp"^,^ in 



the accreting millisecond pulsar IGR J17511-3057 (Falanga 
et al.|[20TT I , and it is indicative of a fuel with a significant 



H fraction, where the proportionality expected is roughly 
tree oc m^^, while the steeper index initially observed in 



T5X2 agrees with models for pure He burning (Gumming 
|fc Bildsten||2000| . In the RB there seems instead to be a 
transition from a fiatter to a steeper regime, but the data is 
too limited to measure the slope of the latter. 

All bursts have a values between 35 ± 5 and 153 ± 26, 
with a weighted average a = 44 ± 14. Values of a for the 
cooling bursts are indicative of a fuel composition that is 
compatible with the solar H abundance, in line with the 
observed long burst durations typical of the H burning rp- 
processes. The higher a values that are measured at the 
highest levels of Fpers are unlikely to be due to a smaller 
H fraction in the burst fuel, as this is the opposite of what 
would be expected for the burst with the shortest tree- In- 
stead, the sudden rise in a might be due to the onset of 
semi-stable He burning in between the bursts. If some frac- 
tion of the fuel is burning in between bursts, this will make 
overall less material available to burn in a burst, therefore 
reducing the fluence and making a grow larger. 

This explanation seems to be backed up by theoret- 
ical modelling of the burst time-scales and energetics in 



this transitional regime. Keek et al. (20121 studied mod- 



els of the NS envelope after a superburst. Bursting is ini- 
tially quenched because of the hot envelope; the first bursts 
reappearing after the superburst are less energetic, feature 
a slower rise and a faster decay than the bursts before it. 
As the envelope cools down further, tree become longer, trisc 
shortens, tdur lengthens and the bursts become again as en- 
ergetic as they were before. We suggest that the same tran- 
sition could be taking place in the RB as the mass accretion 
rate decreases during the outburst decay. Unfortunately, no 
observations are available for the time prior to the appear- 
ance of the non-cooling bursts to look for stable burning or 
the mHz QPOs expected from marginally stable burning at 
higher mass accretion rates. 

It is also worth pointing out that we are not looking at 
the so-called short-waiting-time bursts described by |Keek| 
et al. (20101, which are thought to be due to incomplete 



burning of the fuel in the preceding burst, leading to rapidly 
following, smaller bursts, in sequences of up to four events. 
Within such sequences, the bursts following the first one 
are on average less bright, cooler, and less energetic than 
the initial burst, and their decay profiles lack the longer 
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decay component from the rp-process, suggesting they take 
place in a hydrogen-depleted layer. These short-waiting-time 
bursts show a ~ 5 ( Boirin et al.|[2007 l, a value so low that 
the fuel must be left over from the preceding burst. We there- 
fore exclude this incomplete-burning scenario to explain the 
very short tree we report in the RB, given the stability of 
the burst properties across multiple events, the long burst 
tails indicating the presence of H-burning processes and the 
aforementioned a values typical of ordinary type I bursts 
and the largest in the sample. 



5.3 Bursts at the highest M: indications for a 
slow spin? 

That the bursting rate keeps increasing over such a large 
range in -Fpors is a striking similarity between the RB and 



T5X2. Such behaviour, although predicted by theory (Fuji- 
moto et al.|1981 Heger et al.|200f I has so far been reported 
in no other source. 

When analysing a sample of BeppoSAX Wide Field 
Camera (WFC) data on nine X-ray bursters, both tran- 
sient and persistent sources in the Galactic centre, |Cor-| 
nelisse et al. ( 2003 1 found the burst rate to peak at L ~ 
(1.4-2.1)xl0'''''ergs"^ before dropping by a factor of five. 
Above this threshold, all sources also showed significantly 
shorter bursts, recurring less regularly. It is difficult to in- 



terpret this change in behaviour. Cornelisse et al. (20031 



suggested that it might be due to the transition from un- 
stable to stable H burning, so that bursts above the critical 



luminosity are in a pure Ife regime (the so called case 2, Fu- 
[jimoto et al.|[T98T| , where they take a longer time to build 
up and have much shorter (< 10 s) tails due to the ab- 
sence of prolonged H burning via the rp-processes (LPT93). 
He burning has a much stronger T dependence, meaning 
local conditions set by perturbations in the layer steadily 
burning H have a larger importance in setting ignition con- 
ditions, leading to variations in tree- This seems however an 
unlikely explanation. Firstly, a values should increase sig- 
nificantly while switching to the less efficient burning of He, 
contrary to the result of G08, who reported a values above 
the threshold to stay constant or decrease slightly. Secondly, 
after an initial drop, the bursting rate should start again in- 
creasing through the transition from case 2 to case 1 (the 
He triggered, mixed H/He bursts), again contrary to what is 
observed. Thirdly, case 1 bursts have been observed at lower 
mass accretion rates (e.g. in GS 1826—24, [Galloway et al.| 
2004), suggesting the transition from case 3 to case 2 takes 



place at lower M. Another possibility is that this threshold 
corresponds to the onset of semi-stable He burning (e.g. |in 't| 
Zand et al.|2003 1 , which starts depleting the layer of nuclear 
fuel. The theory of thermonuclear burning however places 
such a transition at near- or above Eddington luminosities, 
inconsistent with the low luminosity (around 10% of L^dd) 
at which it is observed ( Cornelisse et alT| |2003 G08). A fi- 
nal possibility involves variation in burning area. As pointed 



out by Bildsten (20001, the observed flux only reflects the 



overall mass accretion rate M, while the burning regime de- 
pends on the mass accretion rate per unit area m = M /Acb, 
with AbS the surface area of the igniting region. If this area 
increases with a stronger than linear dependence on M, m 
might actually decrease with increasing persistent flux. 

T5X2, by contrast, is very different and behaves more in 



line with theoretical expectations. L12 reported that T5X2 
showed a constant increase in burst rate as its persistent 
luminosity rose from 0.1 to 0.5 I/Edd, with bursts gradually 
evolving into a mHz QPO (with a period as short as 240 s) 
until the onset of stable burning. They posited that this (up 
to then) unique match between theory and observations was 
due to the very slow spin of T5X2. Spinning at only 11 Hz, 
T5X2 rotates an order of magnitude more slowly than all 
other bursters whose rotation rate has been measured via 
burst oscillations. At 11 Hz, rotation is only minimally dy- 
namically relevant, whereas at 100 Hz and above its effects 
are much stronger ( Watts|2012 I. Rotation should, for exam- 
ple, affect the burning process (and its stability) by trigger- 
ing turbulent mixing ( Keek et al.|2009 1. This would stabilize 
He burning at lower accretion rates per unit area than for a 
non-rotating source, skewing the theoretical prediction. 

Given the relatively strong B fleld in T5X2, one might 
wonder whether it is the magnetic fleld rather than the rota- 
tion rate that determines the unusual bursting behaviour in 
this source. Although the fleld is weak enough that fuel con- 
flnement is unlikely to be effective ( Brown & Bildstenll998 1 , 



Cavecchi et al. (20111 argued that dynamical magnetic ef- 



fects may act to conflne the spreading flame front after igni- 
tion, leading to stalling and explaining the presence of burst 
oscillations in this source. Magnetic channelling should also 
lead to a reduction in shearing instabilities and rotational 
mixing in the surface layers independent of the slow rotation 
(L12). Without further study, it has been hard to distin- 
guish these two possibilities for T5X2. However the RB may 
provide a clue, since in this source there is no evidence for 
magnetic channelling (in the form of persistent accretion- 
powered pulsations). 

If rotation rate is the main factor determining the burst- 
ing behaviour in T5X2, then the similarity exhibited by the 
RB implies that it too may rotate slowly. |Patruno et al.| 
( [2012 1 argued that T5X2 has entered its spin-up phase in 
an exceptionally recent time (a few 10^ yr), and that also 
the total age of the binary is extremely low (< 10* yr). It 
therefore has not only a slow spin, but a lso a relatively high 
magnetic field (bet ween 10** and G, |Papitto et al.|2011 



[Miller et al.||201l] ). If the RB also rotates slowly, one may 
speculate that it too is at a similar evolutionary point, a 
young LMXB with a strong magnetic field. Where it differs 
from T5X2, however, is that it also exhibits Type II bursts. 
Interestingly, the only other source known to exhibit type 
II bursts is the slow pulsar GRO J 1744— 28, which spins at 
just 2.14 Hz (Finger et al. 19961. It has a reported mag- 



netic field in the ran ge (2 - 20) x lO"' G ( [Sturner fc Per 
mer[[1996 Cui||1997 l, intermediate between the weak fields 



{B iTTlP - lO^*^ thought to reside in LMXBs and the 
strong fields {B ~ 10^^ - lO" G) found in X-ray pulsars. 
The B value for GRO J 1744- 28 is probably large enough to 
confine burning to a stably burning portion of the surface, 
explaining the lack of type I bursts. 

Confirming the spin rate of the RB is unfortunately dif- 
ficult, since to date there have been no observations of either 
accretion-powered pulsations or burst oscillations from this 
source. A low-significance claim of a 306 Hz modulation in 
a set of stacked burst spectra was reported by [Fox et aT] 
(20011. However results from stacked burst spectra have in 



the past proved misleading, see for example the case of EXO 
0748-676: stacked burst spectra suggested a spin of 45 Hz 
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( Villarreal fc Strohmayer|2004 1 whereas the true burst oscil- 
lation frequency is now known to be 552 Hz ( Galloway et al. 



2010). G08 have also already searched for burst oscillations, 
but only for frequencies above 10 Hz. We are planning a 
more thorough and dedicated investigation, to look into the 
overall timing behaviour of the source, with particular re- 
gard to low-frequency phenomena. (Bagnoli et al., in prep.). 

The lack of persistent pulsations or magnetically- 
confined burst oscillations could be explained by scattering 
in an optically thick Comptonizing cloud, which might de- 
grade the signal and hide the modulation ( [Titarchuk, Cui| 
& Wood 20021, although this seems to be an unlikely ex- 



planation for at least three other LMXBs (G6gii§, Alpar & 
|Gilfanov||2007[ ). Alternatively, no pulsations would be ob- 
served if the dipole moment of the magnetic field were very 



closely aligned with the rotation axis of the NS ( Chen & Ru- 
|derman[T993j . Indeed this would be supported by models for 
type II bursts (e.g. [Spruit fc Taam|1993 l that involve unsta- 
ble disc-magnetosphere interactions in which matter builds 
up at the magnetospheric radius until a critical gas pressure 
is reached. This would not be possible in a highly inclined 
dipole model, where matter can always stream directly to 
the star along a field line. More stringent constraints come 



by D'Angelo & Spruit (2011 20121, who argue episodic 



accretion occurs because of an inner disc radius hovering 
around the corotation radius (the radius at which the 
Keplerian frequency equals the NS spin). They notice that 
the kind of instability they model would not be possible in 
case of a strongly misaligned dipole, where at different lon- 
gitudes the disc would be truncated inside and outside Tc, 



thereby always allowing for some accretion to occur ( Perna, 
IBozzo 



Stella_ 2006) and preventing the buildup of a mass 
reservoir to release later in a type II burst. Although such 
a strongly aligned configuration might seem unlikely, espe- 
cially for relatively young systems, this is in agreement with 
the fact that we only know two type II burst sources, com- 
pared to over a hundred type I bursters. T5X2 could perhaps 
not show type II bursts because it is a more strongly mis- 
aligned rotator than the RB and the bursting pulsar. 
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